Nuclear superfluidity for antimagnetic rotation in 105 Cd and 106 Cd 

Zhen-Hua Zhang, 1 Peng- Wei Zhao, 1 Jie Meng, 1 ' 2,3 Jin-Yan 
Zeng, 1 En-Guang Zhao, 1 ' 4 ' 5 and Shan-Gui Zhou 4,5 

1 State Key Laboratory of Nuclear Physics and Technology, 
School of Physics, Peking University, Beijing 100871, China 
2 School of Physics and Nuclear Energy Engineering, 
Beihang University, Beijing 100191, China 
^Department of Physics, University of Stellenbosch, Stellenbosch 7602, South Africa 
4 State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, 
Chinese Academy of Sciences, Beijing 100190, China 
5 Center of Theoretical Nuclear Physics, 
National Laboratory of Heavy Ion Accelerator, Lanzhou 730000, China 

(Dated: May 17, 2013) 

Abstract 

The effects of nuclear superfluidity on antimagnetic rotation bands in 105 Cd and 106 Cd are in- 
vestigated by the cranked shell model with the pairing correlations and the blocking effects treated 
by a particle-number conserving method. The experimental moments of inertia and the reduced 
B{E2) transition values are excellently reproduced. The nuclear superfluidity is essential to repro- 
duce the experimental moments of inertia. The two-shears-like mechanism for the antimagnetic 
rotation is investigated by examining the shears angle, i.e., the closing of the two proton hole 
angular momenta, and its sensitive dependence on the nuclear superfluidity is revealed. 

PACS numbers: 21.60.-n; 21.10.Re; 23.20-g; 27.60.+j 
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I. INTRODUCTION 



Most rotational bands are usually observed in nuclei with substantial quadrupole defor- 
mations. In these bands, the states decay by strong electric quadrupole (E2) transitions 
and the energy spectra show a pronounced rotational character. Such bands are usually 
interpreted as the result of a coherent collective rotation of many nucleons around an axis 
perpendicular to the symmetry axis. 

In the 1990's, however, a new type of rotational band with strongly enhanced magnetic 
dipole (Ml) transitions and very weak E2 transitions has been discovered experimentally 
in nearly spherical light Pb isotopes [|TJ. This new type of rotational bands have been 
discovered experimentally in a number of nearly spherical nuclei in A = 80, 110, 135, and 
190 mass regions [0411 • The interpretation of such bands in terms of the shears mechanism 
was firstly given in Ref. |J. In order to distinguish this kind of rotation from the usual 
rotation in well-deformed nuclei, the term "magnetic rotation" (MR) was introduced |J. 

In analogy to the antiferromagnetism in condensed matter physics, a similar phenomenon 
known as "antimagnetic rotation" (AMR) is predicted in nuclei by Frauendorf [|[ |7j. The 
AMR band can be explained by the two-shears-like mechanism: in some specific nearly 
spherical nuclei, two valence protons (neutrons) are aligned back to back in opposite direc- 
tions, nearly perpendicular to the orientation of the total spin of the neutrons. A rotational 
band can be built on such a near-spherical nucleus since the rotational symmetry is violated 
by the nucleon currents. Higher angular momenta is obtained by simultaneously aligning 
the two angular momenta for the valence protons (neutrons) toward the neutron (proton) 
angular momentum vector. 

AMR is expected to be observed in the same mass regions as MR ||. However, it differs 
from MR in two aspects. First, there is no Ml transition in the AMR band since the 
transverse magnetic moments of the magnetic subsystems are antialigned. The resulting 
transverse magnetic moment is zero. Second, as the antimagnetic rotor is symmetric with 
respect to a rotation by tt about the total angular momentum axis, the AMR bands consist 
of regular sequences of energy levels differing in spin by 2H and are connected by weak 
E2 transitions reflecting the nearly spherical core. Moreover, the phenomenon of AMR is 
characterized by a decrease of the B(E2) values with spin. Since AMR was proposed 0, it 
has been investigated both from experimental and theoretical aspects. To date, experimental 
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evidence of AMR has been reported in Cd isotopes including 105 Cd @, 106 Cd @, 108 Cd JlO 
110 Cd ||12|| . Most recently, two AMR bands in a single nucleus are firstly observed in 
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107 Cd [13]. In addition, the occurrence of this phenomenon still needs further investigation 
by lifetime measurements in 109 Cd PJ, 100 Pd |5|, 144 Dy [|T§, 101 Pd [T7|], and 112 In 



Theoretically, AMR has been discussed by simple geometry in the classical particle ro 
tor model 



|, and the tilted axis cranking (TAC) model |]19|-pl[]. Based on the TAC 
model, many applications have been carried out in the framework of microscopic- macroscopic 
model P, |T0|, |T5l, and pairing plus quadrupole model P, [T3}|. Most recently, the TAC model 
based on the covariant density functional theory is used to investigate the AMR |22|-f24|[ with 
the point coupling effective interaction PC-PK1 [^5|], for its review see Ref. H. The qual- 
ity of the cranking approximation for principal-axis cranking ||26|| , tilted-axis cranking |27f . 
and aplanar tilted-axis cranking [EH] has been discussed and tested within the particle rotor 
model. 



In the present work, the cranked shell model (CSM) with the pairing correlations treated 
by a particle-number conserving (PNC) method |29], |30] is used to investigate the AMR bands 
in 105 ' 106 Cd. The PNC-CSM method is proposed to treat properly the pairing correlations 
and the blocking effects, and it has been applied successfully for describing the odd-even 



differences in moments of inertia (MOI's) [|31|, the nonadditivity in MOI's [52]], the micro- 
scopic mechanism of identical bands [[$3], the non-existence of nuclear pairing phase 
transition []3"5] , etc. The high-spin states and high-i^ isomers in the rare-earth, the actinide 



region and superheavy nuclei are also well described in the PNC-CSM scheme |36|-[42|]. In 
contrary to the conventional Bardeen-Cooper-Schrieffer (BCS) or Hartree-Fock-Bogolyubov 
(HFB) approach, the Hamiltonian is solved directly in a truncated Fock-space in the PNC 



method [43|. Therefore, the particle-number is conserved and the Pauli blocking effects are 
taken into account exactly. The PNC scheme has been used both in relativistic and non- 
relativistic mean field models [44, 45[] in which the single-particle levels are calculated from 



the self-consistent mean-field potentials instead of the Nilsson potential. 



II. THEORETICAL FRAMEWORK 



The cranked shell model Hamiltonian of an axially symmetric nucleus in the rotating 
frame can be written as 

Hcsm = H + Hp = H N n — ujJx + Hp , (1) 

where H^n is the Nilsson Hamiltonian, —uJ x is the Coriolis interaction with the cranking 
frequency u about the x axis (perpendicular to the nuclear symmetry z axis). Hp is the 
pairing interaction, 

Hp = —G a^a^a^a v , (2) 

where £ (fj) labels the time- reversed state of a Nilsson state £ (77), G is the effective strength 
of monopole pairing interaction. 

Instead of the usual single-particle level truncation in conventional shell-model calcula- 
tions, a cranked many-particle configuration (CMPC) truncation (Fock space truncation) is 
adopted [[50, [46]]. By diagonalizing the Hqsm in a sufficiently large CMPC space, sufficiently 



accurate solutions for low-lying excited eigenstates of Hcsm can be obtained. An eigenstate 
of H CSM can be written as 

|*) = J^Cili) , (a real), (3) 

i 

where \i) is a CMPC (an eigenstate of the one-body operator H ). The expectation value 
of a one-body operator O = X/fcLi is thus written as 

wo\y) = y,chw) +2^2^ (i\o\ 3 ) . (4) 

i i<j 

As O is a one-body operator, the matrix element (i\0\j) for i 7^ j is nonzero only when 
) and \j) differ by one particle occupation [3C|. After a certain permutation of creation 



operators, \i) and \j) can be recast into 

\i) = (-l) M ^.-.) , \j) = {-l) M ^\v ■ ■) , (5) 

where the ellipsis stands for the same particle occupation and (— \) Mi ^ v ) = ±1 ac- 

cording to whether the permutation is even or odd. Therefore, the expectation value of O 



can be separated into the diagonal and the off-diagonal parts 

0(pLv)j , (6) 
^(/x) = </x|<%K, (7) 



= ( l i\0\v)Y^-l) M » +M ''>C i C j , (8) 

i<j 

where = Y2 { \Ci\ 2 Pi^ is the occupation probability of the cranked Nilsson orbital and 
Pj M = 1 (0) if \fi) is occupied (empty) in 

The kinematic moment of inertia of |W) can be written as 

JM = -(*\J X \V) . (9) 

UJ 

The B(E2) transition probabilities can be derived in the semi classical approximation as 

B(E2) = ^|^ol*) 2 , (10) 
where Q20 correspond to the quadrupole moments of protons and 



J - r 2 ) = r 2 Y w . (11) 



III. RESULTS AND DISCUSSIONS 



The Nilsson parameters (k and fi) of 105 > 106 Cd are taken from the Lund systematics [47 . 
The quadurpole deformation parameters are taken as e<i = 0.12 and 62 = 0.14 for 105 Cd 
and 106 Cd, respectively. These values are close to those used in the TAC calculations based 



on the microscopic-macroscopic model or the covariant density functional theory [|9], |22 |. 
The valence single-particle space in this work is constructed in the major shells from N = 
to N = 5 both for protons and neutrons, so there is no effective charge involved in 
the calculation of the B(E2) values. The effective pairing strengths can, in principle, be 
determined by the odd-even differences in nuclear masses and the MOI's, and are connected 
with the dimension of the truncated CMPC space. The dimensions of the CMPC space are 
about 1000 both for protons and neutrons. The corresponding effective pairing strengths 
used in this work are G p = 0.45 MeV and G n = 0.80 MeV for 105 Cd, G p = 0.45 MeV and 
G n = 0.45 MeV for 106 Cd. The data show that the MOI's for the AMR band in 105 Cd are 
smaller than those in 106 Cd. Therefore, a larger effective neutron pairing strength for 105 Cd 



is adopted. A larger CMPC space with renormalized pairing strengths gives essentially the 
same results. In addition, the stability of the PNC-CSM calculation results against the 
change of the dimension of the CMPC space has been investigated in Refs. [RO, [31, [41] . In 



the present calculations, almost all the important CMPC's (with the corresponding weights 
larger than 0.1%) are taken into account, so the solutions to the low-lying excited states are 
accurate enough. 
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FIG. 1. (Color online). The cranked Nilsson levels near the Fermi surface of 106 Cd for (a) 
protons and (b) neutrons. The positive (negative) parity levels are denoted by blue (red) lines. 
The signature a = +1/2 (a = —1/2) levels are denoted by solid (dotted) lines. The deformation 
parameter £2 = 0.14. 

The cranked Nilsson levels near the Fermi surface of 106 Cd for (a) protons and (b) neutrons 
are given in Fig. [j]. The positive (negative) parity levels are denoted by blue (red) lines. 
The signature a = +1/2 (a = —1/2) levels are denoted by solid (dotted) lines. Because the 
Nilsson levels of 105 Cd are very similar with those of the 106 Cd, we do not show them here. 
It can be seen from Fig. [I] that the two proton holes for 105 > 106 Cd are 7r9/2 + [404] (7rg 9 / 2 ). The 
data show that the AMR bands in 105 Cd and 106 Cd are the lowest lying negative and 
positive parity band, respectively. The lowest lying negative parity band for 105 Cd in our 
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calculation is ul/2 [550] (/in/2) an d the lowest lying positive parity band for 106 Cd is the 
yrast band. Therefore, in the following investigation, we do nonadiabatic calculations for 
the 1/1/2- [550] band in 105 Cd and the yrast band in 106 Cd. 




/zoo (MeV) 

FIG. 2. (Color online). The experimental (solid circles) and calculated kinematic MOI's with 
(solid black lines) and without (dashed red lines) pairing correlations for (a) 105 Cd and (b) 106 Cd. 

Figure shows the experimental (solid circles) and calculated kinematic MOI's with 
(solid black lines) and without (dashed red lines) pairing correlations for 105 Cd (upper panel) 
and 106 Cd (lower paner). The pairing interaction is very important in reproducing the 
experimental MOI's, especially the upbending. It can be seen that the MOI's of 105 > 106 Cd 
are overestimated when the pairing interaction is switched off, while they are well reproduced 
after considering the pairing correlations. The first backbending in 105 Cd at fko ~ 0.4 MeV 
is caused by aligning one neutron pair vg-r/2- The configuration after backbending in 105 Cd 
is thus vhxi^ig-j 12) 2 coupled to a pair of 7rg 9 /2 proton holes, which is consistent with the 
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previous calculations || p^j . The first backbending in 106 Cd at fko « 0.3 MeV is caused by 
one pair of neutrons jumping from 1^7/2 to vhn/2- 
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FIG. 3. (Color online). The experimental (solid circles) and calculated B{E2) values with (solid 
black lines) and without (dashed red lines) pairing correlations for (a) 105 Cd and (b) 106 Cd. The 
blue dotted line in (b) is the calculated results with a reduced deformation of £2 = 0.12 in which 
the pairing in considered. The data for 105 > 106 Cd are 

taken from g, g. 



One of the typical features of AMR is the weak E2 transitions reflecting the small de- 
formation of the core. Moreover, the corresponding B{E2) values rapidly decrease with 
the angular momentum, which is connected with the interpretation of the two-shears-like 
mechanism. Figure shows the experimental (solid circles) and calculated B(E2) values 
with (solid black lines) and without (dashed red lines) pairing correlations for 105 Cd (upper 
panel) and 106 Cd (lower panel). It can be seen that the decreasing tendency of the B(E2) 



S 



values with the cranking frequency can be obtained no matter the pairing correlation is 
considered or not. However, the agreement between the data and the calculated results is 
further improved by taking the pairing correlation into account, especially for the higher 
rotational frequency region. For 105 Cd, with paring correlations, the expectation value of 
Q20 decreases from 0.55 eb to 0.41 eb with the rotational frequency huj increasing from 
0.45 MeV to 0.75 MeV. The Q20 value and the corresponding B(E2) value are reduced to 
about 75% and 55%, respectively, which are caused by the effect of the cranking. For 106 Cd, 
it is difficult to describe the B(E2) behavior with a frozen deformation parameter. This 
may be due to the deformation change with the rotational frequency for 106 Cd. In fact, as 
show in the blue dotted line in Fig. |3](b), in order to reproduce the B(E2) behavior from 
huj = 0.45 MeV to hu = 0.75 MeV, a corresponding deformation change from 62 = 0.14 to 
£2 = 0.12 is necessary. Therefore, it can be seen that the two-shears-like mechanism alone 
can provide the decrease of the B(E2) values in 105 Cd, while additional reduction of the 
deformation is needed for 106 Cd. 
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FIG. 4. (Color online). Angular momentum vectors of neutrons J u and the two irgg/2 proton holes 
j n , (a) without pairing (b) with pairing, at rotational frequencies from 0.5 to 0.7 MeV in 105 Cd. 



In order to examine the two-shears-like mechanism for the AMR band, we show the 
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angular momentum vectors of neutrons J v and the two 7rgg/ 2 proton holes j n at rotational 
frequencies from 0.5 to 0.7 MeV in 105 Cd in Fig. It should be mentioned that, in the 
principal axis cranking model, the expectation value of J z vanishes due to the conservation 
of signature. In the present AMR bands, the two proton holes in both 105 Cd and 106 Cd are 
paired. This means that the total angular momentum projection K of these two proton holes 
should always be zero. The angular momenta of the two proton holes could, in principle, 
be extracted exactly from the TAC calculation. Here, we calculate J z approximately in the 



following way according to Ref. [£7 



j z = vwim ■ (12) 

This method has been proved to be a good approximation by comparing the principal axis 
cranking with the particle rotor model in Ref. P7|j . It can be seen from Fig. [|that the two 
proton angular momentum vectors j n are pointing opposite to each other and are nearly 
perpendicular to the vector J v at fiuj = 0.5 MeV. The abrupt increasing of neutron angular 
momentum alignment from hw = 0.5 to 0.6 MeV in Fig. [|(a) is due to level crossing. After 
considering the nuclear superfluidity, the level crossing is delayed and the neutron angular 
momentum alignment increases gradually, which is consistent with the data. With increasing 
cranking frequency the gradual alignment of the vectors j n of the two irgg/2 proton holes 
toward the vector J v generates angular momentum while the direction of the total angular 
momentum stays unchanged. This leads to the closing of the two shears. The two-shears- 
like mechanism can thus be clearly seen. It should be noted that the closing of the two 
proton hole angular momenta becomes more obvious when the pairing correlation is taken 
into account. This indicates the important role played by the the nuclear superfluidity in 
AMR. 



IV. SUMMARY 

In summary, antimagnetic rotation bands in 105 Cd and 106 Cd are investigated by the 
cranked shell model with the pairing correlations treated by a particle-number conserving 
method and the blocking effects taken into account exactly. The experimental moments of 
inertia in 105 Cd and 106 Cd are excellently reproduced with the proper treatment of the nuclear 
superfluidity. The reduced B(E2) transition depends on the deformation rather than the 
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superfluidity. The calculated B(E2) values in 105 Cd are in good agreement with the data. 
In order to reproduce the B(E2) values in 106 Cd, a corresponding deformation change is 
necessary. The two-shears-like mechanism for the antimagnetic rotation is investigated and 
its sensitive dependence on the nuclear superfluidity is revealed. 
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